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Abstract 

We have investigated the stability limits of small spin-polarized clusters consisting of up to 
ten spin-polarized tritium TJ, atoms and the mixtures of TJ. with spin-polarized deuterium DJ, 
and hydrogen HJ. atoms. All of our calculations have been performed using the variational and 
diffusion Monte Carlo methods. For clusters with DJ. atoms, the released node procedure is used 
in cases where the wave function has nodes. In addition to the energy, we have also calculated 
the structure of small clusters using unbiased estimators. Results obtained for pure TJ. clusters 
are in good accordance with previous calculations, confirming that the trimer is the smallest spin- 
polarized tritium cluster. Our results show that mixed TJ.-HJ. clusters having up to ten atoms are 
unstable and that it takes at least three tritium atoms to bind one, two or three DJ. atoms. Among 
all the considered clusters, we have found no other Borromean states except the ground state of 
the TJ. trimer. 

PACS numbers: 67.65. +z, 02.70.Ss, 36.40.-c 
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I. INTRODUCTION 



Electron spin-polarized hydrogen (Hj) has extreme quantum nature due to its very small 
mass and weakly attractive potential. It remains in gas phase even in the limit of zero 
temperature and up to pressures of about 170 bar,- after which it solidifies. Stwalley and 
Nosanow 2 suggested in 1976 the use of H j for achieving a Bose-Einstein condensate (BEC) 
state, which was finally achieved for hydrogen in 1998 by Fried et al.— Heavier isotopes of 
hydrogen, spin-polarized deuterium (Dj) and tritium (Tj) also show a remarkable quantum 
behavior. Dj atoms, which obey Fermi statistics, have nuclear spin one and thus the zero- 
pressure state of bulk D j depends on the number of occupied nuclear spin states. Previous 
theoretical estimations 4 - 1 ^ have shown that (Dji) with only one occupied nuclear spin state 
is gas at zero pressure, while bulk Dj with two (Dj 2 ) and three (Dj 3 ) equally occupied 
nuclear spin states remains liquid at zero pressure and zero temperature. Spin-polarized 
tritium, which obeys Bose statistics, is expected to be liquid-^ 19 - due to its larger mass. 

Recently, Blume et al— have studied microscopic properties of tritium (Tj)jv clusters 
with up to iV=40 atoms using the diffusion Monte Carlo (DMC) method and compared them 
to bosonic 4 HeAr clusters. In that work, it is shown that tritium clusters (Tj)jv are even 
more weakly bound than helium clusters 4 Hejv with the same number of atoms. Furthermore, 
their results suggest the use of TJ. as a new BEC gas with the advantage of a nearly exact 
knowledge of its interatomic potential and the possibility of manipulating the strength of the 
interactions via Feshbach resonances. In addition, it has been shown 1 ^ that the trimer (T j) 3 
is the smallest spin-polarized tritium cluster with a ground-state energy of only —4.2(7) mK. 
It is thus an example of Borromean or halo state,— ^ because the two-body T j system does 
not have a bound state. The prediction of a total zero angular momentum bound state of 
the spin-polarized atomic tritium trimer has been recently confirmed by Salci et al— using 
the finite element method. Importantly, the use of this accurate few-body approach leads to 
conclude that there are not any other bound or shape resonant states with zero and nonzero 
angular momentum.— 

Unlike the case of pure 3 He, 4 He and mixed 3 He- 4 He clusters, which have been the subject 
of both experimental and theoretical study— iI£iI£iILi& nothing is known about the stability 
of mixed or other pure spin-polarized hydrogen clusters different from tritium ones. 3 He 
has almost the same mass as Tj, but because of its small mass and Fermi statistics it takes 
at least 30 3 He atoms to form a bound state.— Recently, Hanna and Blume^ 9 - have studied 
the near-threshold behavior of weekly bound three-dimensional bosonic clusters with up to 
40 atoms interacting additively through two-body van der Waals potentials modeled after 
the Tj-Tj interaction. Among other results, they found a four-parameter fit for the critical 
mass needed to bind a bosonic cluster of N atoms. From this, it follows that if Dj were a 
boson, due to its smaller mass compared to Tj., it would take at least 13 atoms to form a 
bound cluster (Dj)^. Considering that Dj obeys Fermi statistics we expect that, similarly 
to the case of 3 He, 18 it will take more than 30 Dj atoms to form a self-bound cluster. On 
the other hand, although Hj obeys Bose statistics, its approximately three times larger 
zero-point energy (with respect to Tj) implies that clusters of Hj are not stable. 19 This 
statement is in agreement with the gas nature of bulk Hj at zero temperature.- Therefore, 
hydrogen-tritium and deuterium-tritium clusters are reasonable choices for getting stable 
self-bound spin-polarized clusters. 

If tritium could be experimentally prepared to achieve a BEC state, as suggested in 



Ref. llOj, the presence of a finite fraction of deuterium would be a new example of a Bose- 
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Fermi mixture with BEC. In fact, Bose-Fermi mixtures composed with alkalines are nowa- 
days studied in several laboratories and the main physics underlying them is well known.— 
Tritium-deuterium mixtures would be even more interesting since they are isotopic and with 
a very well known interaction. Dilute vapors are also candidates for searching the elusive 
Efimov states,— and recently experimental evidence of its existence has been observed in an 
ultracold gas of cesium atoms.— In the case of pure tritium, it has been proved that such 
state does not exists but mixed isotopic hydrogen systems could offer new possibilities for 
continuing the search of trimers with Efimov character. 

In this article, we report energetic and structural properties of small spin-polarized 
tritium-deuterium and tritium-hydrogen clusters, obtained using the diffusion Monte Carlo 
method. It is our main goal to determine which of these clusters are stable. As many of 
them extend beyond the classically forbidden region they can be considered as quantum 
halo systems. We are particularly interested in finding among these halo clusters mixed 
Borromean or " super-Borromean"— systems for which all subsystems are unbound. So far, 
3 He2K is the only mixed molecular system for which the Borromean state has been predicted 
theoretically.— 

In Sec. II, we briefly describe the DMC method and the trial wave functions used for 
importance sampling. Due to the fermionic nature of DJ, , the fixed-node and released-node 
methods are used in cases where the wave function has nodes. Sec. Ill reports the results 
obtained, and finally, Sec. IV comprises the summary of the results and main conclusions. 



II. METHOD 

We consider that all the HJ,, DJ,, and TJ, atoms interact with the spin-independent central 
triplet pair potential & 3 £„. This interaction has been determined in an essentially exact 
way by Kolos and Wolniewicz,— and recently extended to larger interparticle distances by 
Jamieson et al. (JDW).^ We have used a cubic spline interpolation between JDW data 
and smoothly connected the points to the long-range behavior as calculated by Yan et al..— 
The same form has been used in the recent DMC calculation of bulk HJ,, 1 work in which 
a comparison between different potentials employed in the past is also reported. On the 
other hand, we have checked that the addition of mass- dependent adiabatic corrections (as 
calculated by Kolos and Rychlewski 27 ) to the JDW potential does not modify the energy of 
spin-polarized hydrogen clusters. 

The starting point of the DMC method is the Schrodinger equation written in imaginary 
time, 

-h^^- = (H-EMR,t), (1) 

where E is a constant acting as a reference energy, R = (ri, . . . , r^) is a walker in Monte 
Carlo terminology, and H is the iV-particle Hamiltonian 

"=-£^V? + I>M. (2) 

i=l 1 i<j 

As usual in the method, a trial wave function ip(R) is introduced for importance sampling 
and then the Schrodinger equation ([1]) is rewritten in terms of <&(R, t) = ^(R, t)ip(R). The 
resulting equation is solved stochastically by considering a short-time approximation for the 
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Green's function. In the limit t — > oo (for long simulation times), only the lowest energy 
eigenfunction not orthogonal to ^(R) survives and then the sampling of the ground state is 
effectively achieved. 

The trial wave function used for the simulation of the pure T j clusters is of Jastrow form, 
i/jj(R) = Yli<j /r( r ij), where describes correlations between pairs of tritium atoms. 

Similarly, wave functions for the study of mixed tritium-hydrogen clusters, containing Ni 
Hj atoms and N 2 Tj atoms (Ni + N 2 = N), are constructed as a product of two-body 
correlation functions between all the pairs, 

ATi N ATi N 

^(K;N 1 ,N 2 )=l[f H (r ij ) J] / T (r„) II fMnj) , (3) 

i,j=l i,j=Ni+l *=1 i=iVi+l 

i<j i<j 

where /n(r) describes the two-body correlations between HJ, atoms and furir) accounts 
for the Hj-Tj pairs. In case of mixed (Dj^^Tj)^ cms t ers > the total wave-function is 
constructed as a product of an antisymmetric function corresponding to deuterium (i^a) 
and Jastrow factors, 

Ni N Ni N 

^(R;iv 1 ,iv 2 )=^n^(^) n n w^o, (4) 

i,j=l i,j=Ni+l »=1 j=Ni+l 

i<j i<j 

where /d(^) describes two-body correlations between Dj atoms and fori?) two-body corre- 
lations for the D j-Tj pairs. For the function ipA it is enough to consider the product of Slater 
determinants for each one of the nuclear spin states. Consequently, ipA has different form in 
the case of Dj 2 and Dj 3 mixed Dj-Tj clusters. For example, in the case of (Dj^^Tj)^ 
clusters, ipA — 1 f° r ^1 — 1,2, ipA= x i — %i for N\ = 3, and ipA=( r i — r 2)( r 3 — r 4) for 
Ni=A, while for (Dl 3 ) Nl (T l) N2 clusters ipA — 1 fo r ^i — 1,2,3, if) a — x\ — x 2 for Ni=4 
and ipA=( r i - r 2 )(r 3 - r 4 ) for N x =5. 

We have worked with three different models for the two-body correlation function /t(^), 

f T (r) = exp[-6i exp(-6 2 r) - b 3 r] , (5) 



1 



r 



f T (r) = exp 



f T (r) = - exp - - - sr , (6) 



— sr 



(7) 



where b\, 6 2 , ^3, 7, s and b are variational parameters. They have been obtained by 
optimizing the variational energy calculated with the VMC method. Finally, we found the 
form ([6]) to be optimal for clusters from three to five T j atoms, while the other two (I51P7I) suit 
better for larger clusters. The same function, with different variational parameters, has then 
been used for all the other two-body correlation functions (/#(?"), /i)(r), /^r(^), /dtW)- 
For example, in case of pure Tj clusters (a, 7, s) range from (3.9,4.0,0.0012) for trimer to 
(4.0,3.7,0.001) for (T|) 5 and (b,s) range from (3.55,0.065) for (Tj) 6 to (3.6,0.05) for (Tj) 10 . 
The VMC recovers from 70% to 90% of the DMC energy, except in case of (Tj) 3 where 
the VMC energy is only -0.85(8) mK. In case of mixed clusters the parameters in the fr{ r ) 
function do not change significantly, so for example in case of the form ([7]) b assumes values 
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from 3.5 to 3.62 and s from 0.07 to 0.05, going from smaller to larger clusters, respectively. 
At the same time, in the function /^(r), b goes from 3.58 to 3.7 and s is around 0.003. 
Similar behavior of variational parameter b is obtained for the other correlation functions, 
while s is around 0.01 in forif) and around 0.003 in /o(r). In most cases of mixed clusters 
VMC recovers from 60% to 90% of the DMC energy, exceptions being smaller clusters with 
more than two DJ, atoms and the clusters having only three Tj atoms, where the VMC 
energies are further from the DMC ones. 

The DMC method we have used is accurate to second order in the time step At.— The 
ground-state energies have been calculated for several time-steps and then extrapolated to 
zero At to remove any possible time-step bias. We have also studied the optimal mean 
walker population and finally chosen N w = 1000 in order to eliminate any bias coming from 
it. 

Apart from statistical uncertainties, the energy of the bosonic clusters is exactly cal- 
culated. In the case of mixed clusters with more than three D| 3 or two D| 2 atoms, the 
sign problem appears because the trial wave function changes sign. Since the Monte Carlo 
method requires that /(R, t) = ^(R, t)^(R) > 0, we have first used the so-called fixed-node 
approximation which allows only the moves in which if) and ^ change sign together, thus 
fixing the nodes. In this way, an upper bound to the energy is obtained.— In a subsequent 
step, we have removed the nodal constraint imposed by the trial wave function by using the 
released-node method. We have adopted the methodology used previously in liquid helium 
calculations,— and recently employed in the study of small mixed helium clusters. 17 In this 
approach, walkers are allowed to cross the nodes imposed by the trial wave function and 
survive for a finite lifetime. In order to achieve an effective crossing of the nodal surface 
an auxiliary guiding wave function if) g is introduced in the DMC calculation. This function 
is positive everywhere, different from zero in the nodes, and approaches \ip\ away from the 
nodes. Like in Refs. Il7|]30l . we have taken ip g = , 0j(R)(^(R) + a 2 ) 1 / 2 , which satisfies the 
above condition for suitable choices of a. The released-node energy is estimated through an 
exponential fit E{t) = E r + Ae~^l T ^ to the DMC data, with t the released time. In all cases 
where the RN method is used in this work the difference between the last calculated point in 
released time and E r is of the same order as the statistical noise. The behavior of a and the 
maximum released time needed to achieve the asymptotic limit of the released energy has 



been similar to the case studied in Ref. [17J. Also, like in Ref. [17|, we find the RN energies 



to be very close to the FN ones. In addition, in the case of (T \)a{D | 2 )3 we have verified 
that the inclusion of backflow correlations in ipA does not change the final FN results. 



III. RESULTS 

DMC results for the ground state energy of pure and mixed Tj-Hj clusters are given 
in Table [I] Our results for pure tritium T j clusters are compared to the results of Blume 
et. al— in Fig. [TJ We obtain a good agreement with these published data, although 
our results are slightly lower for all N. Underlying this difference is the fact that in Ref. 



lOl a damped three-body Axilrod- Teller potential term 31 is introduced in the Hamiltonian, 
raising the trimer energy 1.6% and the energy of the cluster with 40 particles 6%. We also 
find that the trimer is the smallest spin-polarized tritium cluster with an energy of only 
—4.8(2) niK, confirming it to be a halo state. In order to test the sensitivity of the present 
results to the details in the interaction potential, we have also calculated the energy with the 
potential that Silvera and Goldman 3 ^ constructed as a fit to the older Kolos and Wolniewicz 
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data. 24 It is worth mentioning that a full comparison between the different potentials used 
in the literature for studying bulk Hj can be found in the Ref. H|. Using the Silvera model, 
we obtain a trimer binding energy significantly smaller, —1.9(4) mK. The trimer energy 
is therefore sizably affected by the interatomic potential and the reason for that lies on 
the huge cancellation between the kinetic and potential energies. For instance, using the 
JDW potential and the pure estimator for the potential energy, we obtain E p = —355(9) 
mK and E^ = 350(9) mK. For larger clusters, the influence of the potential form on the 
binding strength is reduced: for the tetramer and with the Silvera potential, we obtain 
an energy E = —107(2) mK which is only about 15% weaker than the one for the JDW 
potential, E = —126(2) mK. The potential and kinetic energies of the tetramer using the 
JDW potential are E p = -1871(11) mK and E k = 1745(12) mK. 

The addition of one Hj atom to the (Tj)^ cluster creates a system that appears to be at 
the threshold of stability or unstable, that is the ground-state energy of the (Tj)]vHj cluster 
is within the errorbars equal to or higher than the energy of the (Tj)^ cluster, respectively. 
In order to be certain about the conclusion concerning the stability of the (Tj)jvHj clusters, 
we have repeated the calculations using two different types of trial wave functions. In the 
case where VMC energies are further away from the DMC ones, a lot of simulation time 
is needed for the system to reach the equilibrium state, and in some cases with the trial 
function ([6]) the lowering of the energy with imaginary time going on is almost imperceptible. 
On the other hand, with both ([5]) and ([7]), and for small values of the parameter s or b% 
in the f{rTH) function, the system rather quickly reaches the equilibrium energy, which is 
within the errorbars equal to the energy of the system without the Hj atom. Furthermore, 
if we continue the simulation after the system has reached the ground-state energy, we can 
observe the surplus Hj atom leaving the cluster. This indicates that these clusters are 
effectively unstable. We have also calculated the ground-state energies of several clusters 
with more hydrogen atoms (Tj)7v(Hj)M, such that N + M < 10. Our results indicate 
that these clusters are also unstable, because although in some cases negative energies are 
obtained they are above those for (Tj)jvHj. 

Table HT1 presents the results for the ground-state energy of mixed Tj-Dj clusters. Due 
to the approximately twice larger mass of Dj compared to Hj, clusters with one or two 
Dj atoms and at least three Tj atoms are stable. In the clusters (Tj)jv(Dj) 2 , we have 
assumed that Dj atoms occupy two different nuclear spin states. For clusters (Tj)jy(Dj) M 
with M >2, we have considered both the case where two different nuclear spin states are 
occupied (Tj)jv(Dj 2 )M and the case with three different occupied nuclear spin states (Tj 
)iv(Djs)M. Similar to the case of the bulk system,-^ we find that the three-component 
spin clusters are more strongly bound than the ones where Dj occupies two nuclear spin 
states, £ , ((Tj) A r(Dj 3 )A/) < E((Tl) N (B j 2 )iw")- Furthermore, (Tj)7v(Dj 3 ) 3 , whose spatial 
wave function has no nodes, are all stable for N > 3, while it takes 6 Tj atoms to form a 
stable (Tj)7v(Dj 2 )3 cluster. The smaller (Tj)7v(Dj 2 ) 3 clusters have the same energy as the 
(Tj)jv(Dj) 2 clusters within the errorbars, but the analysis of the distributions shows one of 
the Dj atoms leaving the system. The clusters (Tj)7v(Dj 3 ) 4 are similarly at the threshold 
of binding for N = 3, 4, 5, 6, just like the clusters (Tj)Ar(Dj 2 ) 4 , and again it is necessary to 
consider the structure to better determine their stability. Our analysis indicates that only 
(Tj) 6 (Dj 3 ) 4 is stable. The same behavior, as far as the energy is concerned, is reproduced for 
(Tj)Ar(Dj 3 ) 5 with N < 5: as the separation of the Dj atoms is growing along the simulation 
we are led to consider them unstable. We expect that all other clusters (Tj)Ar(Dj)M, with 
iV + M < 10 and M > 5 are unstable because the exchange of T j with the D j atoms does 
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not change the interaction but raises the kinetic energy. The evolution of the energy of 
mixed clusters formed by N tritium atoms and one hydrogen atom or up to four deuterium 
atoms is shown in Fig. [2l 

In addition to the energy, DMC simulations allow also for exact estimations of other rel- 
evant magnitudes such as the distribution of interparticle distances P{r) or the distribution 
of particles with respect to the center of mass of the cluster p(r) . In both cases it is possible 
to eliminate the bias coming from the trial wave function by using pure estimators^ 3 - and 
arrive to exact results. 

Our calculations confirm that the pure Tj clusters are spatially very diffuse.— From the 
DMC results for P(r) one can see that the average separation between particles (vtt) ranges 
between 34 A in the case of (Tj) 3 to 10.9 A for (Tj) 6 and 10.5 A for (Tj)i . The addition 
of one Dj or one Hj reduces the average separation between Tj atoms in the case of clusters 
with N = 3 or 4, as can be seen in Fig. [31 which compares the pair distribution functions 
of (Tj) 3 , (Tj) 3 Dj and (Tj) 3 Hj. In contrast, for clusters with more Tj atoms the Tj-Tj 
separation remains almost unaffected by the addition of one Dj or one Hj atom. Due to 
the larger zero-point motion of D j and Hj atoms, (r TT ) > (r TD ) > (r TH ). Moreover, we 
obtain (tth) greater than 100 A indicating that the cluster (Tj) 3 Hj is unstable. Similar 
behavior is noticed for larger (Tj)ArHj clusters. 

The effect of adding Dj atoms to the core of (Tj) 6 cluster is shown in Fig. HI The 
separation between pairs of Tj atoms remains almost unchanged, while the Tj-D j and D j- 
Dj separations grow with the addition of D j atoms. There is a noticeable difference between 
clusters with Dj 2 and Dj 3 spin-polarized deuterium. Former clusters with more than three 
particles are much more extended, i.e., the average separation between Dj atoms is larger 
and the tail of the distribution decays slowly. Despite the large size of the cluster, the 
average separation between particles does not appear to grow in the case of (Tj) 6 (Dj 2 ) 3 so 
we are led to consider this cluster as stable; the same conclusion applies to (Tj) 6 (Dj 3 ) 4 . On 
the other hand, the average Dj-Dj distance in the (Tj) 6 (D j 2 )4 cluster is slightly above 40 A 
but grows very slowly in the course of the simulation, indicating that it could be unstable. 

Finally, Fig. [5]presents the density distribution of the clusters (Tj)Ar(Dj), for iV =3,6,9. 
With the increase in the number of bosons, deuterium is pushed to the surface of the 
cluster, although even for N = 9 there is still an appreciable probability of finding it inside 
the cluster. Similar behavior is found for (Tj)jy(Hj) clusters. For larger clusters, one can 
expect that both Dj and Hj will be pushed to the surface forming the so-called Andreev 
states,— similarly to the well known behavior of 3 He in mixed 3 He- He clusters. 15 



IV. CONCLUSIONS 

The ground-state properties of spin-polarized pure and mixed tritium clusters have been 
accurately determined using the DMC method. Our results show that the trimer is a Bor- 



romean or halo state in agreement with the results of Ref. [10j. The most promising can- 
didates for the super-Borromean states in mixed clusters have been (Tj) 2 (Dj)i i2 , 3 ,4,5, but 
after a careful analysis we have found no bound states in either of these clusters. On the 
other hand, we conclude that it takes at least three T j atoms in order to bind 1-3 D j atoms, 
the case of (Tj) 3 (Dj) 3 being bound only when the deuterium atoms are of Dj 3 type. The 
DMC results show that 6 T j atoms are needed in order to bind 4 D j 3 atoms, while clusters 
with up to total 10 atoms and more than 4 Dj atoms seem all to be unstable. We have not 
considered mixed clusters with the D j atoms in only one nuclear spin-state Djx since, from 
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analogy with the bulk, we expect them to be less bound than D| 2 and D| 3 . 

Finally, and concerning mixed tritium-hydrogen clusters, our results show that all the 
(TJ,)jv(HJ,)m clusters with N + M < 11 are unstable. Work is in progress to find out how 
many Tj atoms are needed to bind one Hj atom. 
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Figure captions 

FIG. 1: Comparison of our results for spin-polarized tritium clusters (circles) with the 
results of Blume et. alM (crosses). The difference between our results and the ones from 
Blume et. al. 10 is mainly due to the use of a slightly different interaction (see text). The 
error bars of the DMC energies are smaller than the size of the symbols. 

FIG. 2: Energies of the spin-polarized pure and mixed tritium clusters in K as a function 
of the number of T j atoms. 

FIG. 3: Distribution of interparticle distances in the (Tj) 3 , (Tj) 3 D j and (Tj) 3 Hj clus- 
ters. 

FIG. 4: Distribution of interparticle distances of the (TJ,)e(DJ,)jv clusters. Full line 
denotes Tj-Tj, dashed line Tj-Dj, and dotted line Dj-Dj separation. N = 3,4 distributions 
with lower maxima in Dj-Dj separations correspond to D I2 and those with higher maxima 
toDj 3 . 

FIG. 5: Density distributions of Tj (upper curves) and Dj (lower curves) for (Tj^Dj 
clusters. 
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TABLE I: The ground state energy (in mK) of pure spin-polarized tritium clusters with N atoms 
and tritium clusters with an additional hydrogen atom. Clusters which appear to be unstable are 
written in italic. Figures in parenthesis are the statistical errors. 



E[(Ti) JV (Di) Af ] 
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-2237(6) 
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TABLE II: The ground-state energy E (in mK) of mixed spin-polarized tritium-deuterium clusters. 
For mixed clusters having more than two deuterium atoms, combinations with Dj 2 and D| 3 are 
considered. The clusters who are at the threshold of binding or appear to be unstable are written 
in italic. Figures in parenthesis are the statistical errors. 

1 1 1 1 1 1 1 1 

e 

© 

* 

* 

1 ' 1 1 . . . . 

2 4 6 8 10 

N 

FIG. 1: I. Beslic et al. 
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FIG. 4: I. Beslic et al. 



11 




12 



